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Abstract. 

We have characterized numerically, using the Janus computer, the Lee- Yang 

complex singularities related to the overlap in the 3-D Ising spin glass with binary 
couplings in a wide range of temperatures (both in the critical and in the spin-glass 
phase). Studying the behavior of the zeros at the critical point, we have obtained an 
accurate measurement of the anomalous dimension in very good agreement with the 
values quoted in the literature. In addition, by studying the density of the zeros we 
have been able to characterize the phase transition and to investigate the Edwards- 
Anderson order parameter in the spin-glass phase, finding agreement with the values 
obtained using more conventional techniques. 



PACS numbers: 75.10.Nr,71.55.Jv,05.70.Fh 
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1. Introduction 

In two seminal papers, T. D. Lee and C. N. Yang [H |2] introduced a new tool to 
understand the origin of a phase transition by studying the complex singularities of 
the free energy, or, equivalently, the zeros of the partition function. In particular, they 
showed that all the zeros are located on the unit circumference on the complex activity 
plane (taking as variable z = e~^'^, where h is the external magnetic field). They also 
proved that if the zeros condense onto the real axis when V ^ oo a phase transition takes 
place. Finally, they related, in the low-temperature phase, the density of zeros with the 
discontinuity in the order parameter (remember that the Ising model experiences a first- 
order phase transition when h changes at a fixed temperature below the critical one). 
This approach was subsequently extended to the temperature zeros by Fisher [21 SI El E] • 
We have Lee- Yang like theorems only for a limited class of non-disordered systems 
(such as Ising models). However, it is possible to develop a scaling theory by assuming 
that asymptotically the complex singularities (wherever they lie) touch the real axis 
(thus generating the phase transition). Hence, despite the lack of formal theorems it 
is still possible to apply Lee and Yang's main results to a wide class of systems (e.g.. 
Potts models [7|). In this class of systems the zeros do not live on a circle as stated by 
the Lee- Yang theorem, but they still control the critical properties of the model. We 
will only assume this last fact irrespectively of the form of the locus of the zeros in the 
complex plane [7]. 

In Refs. [HI m [ID] the analysis of complex singularities was applied successfully to 
diluted systems (in particular diluted Ising models in two and four dimensions). The key 
point for the applicability of the standard results, well tested in non-disordered systems, 
is to compute the complex singularities individually for each disorder realization (called 
sample) and then use the mean of the individual zeros (sample zeros) in order to test 
the scaling properties of the zeros and to study the properties of the integrated density 
of zeros. In this work we will also introduce the analysis of the median. 

Nowadays we are interested in gaining a deeper understanding (from the point of 
view of the complex singularities) of the properties of an interesting frustrated and 
disordered system: the three-dimensional Ising spin glass. The magnetization, while 
very interesting in off-equilibrium dynamics and in experiments, plays no role in the 
critical behavior and in the understanding of the low-temperature properties in a finite- 
dimensional spin glass. The observable that controls this spin-glass phase is the overlap. 
Hence, in this work we have focused on the numerical study of the complex singularities 
linked with the overlap in order to study the phase transition and the properties of the 
spin-glass phase. 

In the past, Lee- Yang and Fisher zeros were obtained in spin glasses by means of 
the numerical evaluation of the partition function on small lattices [HI [121 lEl [E] • This 
methodology was also applied to models defined on Bethe lattices [15]. Finally, some 
calculations were performed with the help of replicas [IB] . 

More recently, the complex singularities linked with the external magnetic field 
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were studied for the two and three-dimensional Ising spin glass model in the interesting 
reference |T7], which focuses on the Griffiths singularity and computes all the zeros for 
small lattices. 

In particular we are interested in characterizing the scaling of the individual zeros at 
the critical point (which will allow us to compute the anomalous dimension exponent) 
and checking the scaling in the spin-glass region. In addition, we want to study the 
properties of the density of zeros in the critical and spin-glass region: the behavior 
of this observable will clearly signal the phase transition. Finally, we will show how 
this density of zeros can be used to compute the Edwards-Anderson order parameter. 
However, the spin-glass susceptibility presents strong scaling corrections (even on an 
L = 32 lattice and /3 = 1.4), which induce strong corrections on the density of zeros 
allowing us (from the numerical point of view) only to test our density of zeros against 
the values of gsA found in the literature, rather than attempting a direct numerical 
computation of the order parameter. We want to stress that in cases in which the spin 
glass susceptibility reaches the asymptotic value, the method we propose will be able to 
provide directly the order parameter (^ea) giving an additional method to those used 
nowadays [ISl US EDI EI] ■ 

Let us finally mention that we have obtained the data presented in this work from 
the analysis of the configurations produced in parallel tempering runs [201 EI] using the 
Janus computer [22l|23|[2i]. 

2. Model and observables 

We have studied the three-dimensional Edwards-Anderson model with dynamical 
variables (jj. These variables are Ising spins and are placed on the nodes of a cubic 
lattice of linear dimension L and volume V = L^. The Hamiltonian of the system is 



where {ij) indicates that the sum is over the nearest neighbouring sites. The couplings 
Jij are random quenched constants with bimodal probability distribution, that is, 
J = ±1 with 50% probability. Every realization of the couplings is called a sample. Due 
to the fact that we have a random Hamiltonian, we have to deal with a double average: 
first the thermal average, which we will denote by {{■■■)), and then the average over the 
samples, which we will denote by (•■■). 

We have simulated several real replicas of the system, so we can define the local 
overlap 

«. = ' (2) 

where ai^^ belongs to the first replica and cr-x^ belongs to the second one. The spin 
overlap is defined from this local overlap as 




(1) 




(3) 
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where the sum runs over the whole volume iV). In addition, we define q = Q/V. These 
observables allow us to define some new quantities, for example the non-connected spin- 
glass susceptibility 

1 



XSG=y{Q'). (4) 

Let us now rewrite the Hamiltonian adding a new perturbation eQ and including 
the two replicas explicitly, 

n, = n{a^)+nia2) + eQ 
This Hamiltonian looks like that of the Ising model in a magnetic field 

nh = no + hM. (6) 

We can write the partition function, whose zeros we want to study, as 
= (cos(/3eQ)e-^^o + zsin(/3eQ)e-^^°) . 

[(7(1)<t(2)] 

Let Zq the partition function of the non-perturbed system, so 

Z = Zo{(cos(/3eQ)) + z(sin(/3eQ)} (8) 

and we have to find the zeros of the function (cos(/3eQ)) since in absence of a magnetic 
field {sm{(3eQ)) is zero. The algorithm to find them is quite easy: we start from the 
list of individual measurements of Q for each sample (see Section |4]) and evaluate the 
average (cos(/3eQ)), increasing e in small steps Ae. When the function changes signs 
from one step to the next, we have found a zero in this interval. Obviously, the smaller 
Ae the better the precision of the zero that we have found, but also the slower the 
analysis, so we have to be careful with the error estimates. We have analyzed the first 
four zeros of this function. 

3. Finite-Size Scaling 

One can obtain the expected behavior of the LY zeros by means of (see for example Ref. 

m) 

VxsgV 

therefore, the finite-size dependence, at the critical point, of the Yang-Lee zeros can be 
expressed as: 

e,(L)~L--S (10) 

where 

x, = {D + 2-r])/2, (11) 
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and D is the dimensionality of the system, being D = 3 in this work. If corrections to 
scahng are taken into account, the above relation becomes 

e,{L)^L'^^[l + 0{L--^)) , (12) 

where X2 is the leading correction-to-scaling exponent, X2 = oj. 

In the broken symmetry phase, where the non-linear susceptibility diverges as the 
volume of the system, we expect the following behavior: 

^^L)^^, (13) 

We can take scaling corrections into account, as in the critical point, and 

ej{L) ^ V-^ {I + 0{L-^')) . (14) 

where is the leading correction-to-scaling exponent in the broken phase ||j 

In order to discuss the density of zeros we need to describe some known properties 
of the Hamiltonian defined in Eq. ([s]). This Hamiltonian was introduced in the 
past [23 En]- In particular it experiences a first-order phase transition in e, below 
the critical temperature of the uncoupled model. Hence, the overlap is discontinuous: 

lim (g)(e) = ±gEA, (15) 

being the discontinuity at the transition just 2gEA- 
We can also introduce the density of zeros 

/x,(e) = ^5(e-e,(L)) (16) 

and its integrated version 

G{e) = [ dxfx,{x) (17) 
Jo 

which takes the following value computed for a given zero: 

2 j- l 
2V 

where j labels the j-th zero (j = 1,2,...). In order to deal with the discontinuous 
behavior of G{e) at the zeros, we follow the recipe of references [27[ [2H] and use the 
mean between two consecutive plateau values (j — 1 and j). Anyhow, the asymptotic 
value of the integrated density computed in the j-th zero is j/V. We will discuss this 
point again in subsection |5.3[ 

This integrated density is very useful to characterize a phase transition. In general 
it behaves as 

G(e) = + ag (19) 

I Both droplet and RSB predict algebraic decays for the connected correlation functions in the spin 
glass phase (the spin glass phase is critical in both models). In the droplet model the exponent of the 
decay is y (sometimes denoted as 0), so one can show that — y. In RSB depending of the value of 
q we have different decays (of the g-constrained correlation functions), denoting the decay exponents 
as 9{q). So the leading correction exponent can be shown to be the smallest of the different 0{q). See 
Refs. [iniHI] for a detailed discussion on 9{q). 



G{e,{L)) = ^, (18) 
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and we can extract a great amount of physical information from these three numbers 
(ai, 02 and 03): 

• In the symmetric phase 03 < 0. In a broken phase 03 > 0. 

• In the onset of a first-order phase transition, varying e as it is our case: 02 = 1 
and 03 = 0. In addition we can extract the order parameter of the broken phase: 
gEA = 7rai//3||| 

• At the critical point, 03 = and 02 is related with the anomalous dimension 77 by 
means: 

a2 = — . (21) 

^ + 2-7] ^ ' 

4. Simulation details 

We have run simulations for several lattice sizes on the Janus supercomputer [221 ESI El] 
(for L = 16,24, and 32) and on conventional computers (for L = 8 and 12). These 
simulations were originally reported in [20], which gives full details on the chosen 
parameters and the thermalization protocol. In this section we give only a brief 
summary. 

We have used the parallel tempering algorithm [29], |30] , choosing the temperatures 
to maintain an acceptance around 20% in parallel tempering updates. Besides, since 
Janus needs far more time to do a parallel tempering update than a heat-bath one, we 
have chosen to do one parallel tempering update every 10 heat-bath ones. In table 
[1] one can find a summary of the simulations parameters. In order to choose the 
simulation length, we have assessed thermalization on a sample-by-sample basis, using 
the temperature random walk technique [SD ED] (table [l] gives the average number of 
lattice updates for each L). 

In general, each of the single processors (FPGAs) of Janus takes care of the 
simulation of one replica of the system. However, some samples have such a slow 
dynamics that even with this algorithm the simulation would be too long (more than six 
months of continuous running time), so we would need to accelerate it. For these few 
cases we have created a special low-level code that is in charge of the parallel tempering 
in the control FPGA of a board of Janus. This allows us to spread the simulation 

§ In Lee and Yang's paper, the starting point is the Hamiltonian fSHh = PV-h^o + hM, where M is the 
total magnetization of the system. In terms of the fugacity z = e"^'', they obtained the following result 
(valid below the critical temperature) for the density of zeros (in the fugacity variable that we will 
denote as /i^(z)): /i^(0) = r7isp/(27r), where TOgp is the spontaneous magnetization below the critical 
temperature. In order to transfer this result to our notation we remark that our "magnetic field" is /3e, 
(7ea plays the role of mgp and we need to use the standard law of the transformation of the probability 
densities (/^^(z) = ii^{t)\del dz\, where z — exp(— 2/3e)), obtaining: 

/i.(0) = ^. (20) 

TT 

Notice that near e = we can identify ai with /ic(0). 
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over several processors running only a subset of temperatures in each FPGA, thus 
accelerating the simulation by increasing the parallelism. 

Table 1. Summary of the simulations. A^t is the number of simulated temperatures 
(evenly spaced between Tmin and Tmax); -^mes is the number of Monte Carlo steps 
(updates of the whole lattice) between measurements; is the average simulation 

time (since we use the random-walk technique the simulation time depends on the 
sample); A^sam is the number of simulated samples. We have simulated four real 
replicas for each sample. Finally, L = 8 and L = 12 have been simulated on PCs 
and L = 16, i = 24 and L = 32 on Janus. 



L 


T ■ 

-'- mm 


T 

-'■ max 




N 

^ ' mes 




N 

^ 'sam 


8 


0.150 


1.575 


10 


10^ 


7.82 X 10^ 


4000 


12 


0.414 


1.575 


12 


5 X 10^ 


3.13 X 10^ 


4000 


16 


0.479 


1.575 


16 


10^ 


9.71 X 10^ 


4000 


24 


0.625 


1.600 


28 


10^ 


4.02 X 10*^ 


4000 


32 


0.703 


1.549 


34 


2 X 10^ 


1.90 X 10^° 


1000 



4.1. Data for the computation of the zeros 

We have saved on disk every individual measurement of the overlap. Since we have 
simulated four real replicas of the system, for each sample we have a total of 6A'^HB/A^mes 
values of Q. Given the variable iVnB, this ranges from 1.2 x 10^ to 2 x 10^ measurements 
for our largest lattices, so we have a very good precision for computing the zeros 
of the partition function. We have discarded the first half of the measurements for 
equilibration. 

We want to study the behavior of the system in the critical temperature and in the 
low-temperature phase of the system, analyzing the scaling of the zeros. Therefore, 
we need to compute the zeros of different system linear sizes, L. but at the same 
temperature. Since wc have not simulated the same temperatures for every lattice 
size, we have interpolated, using cubic splines, in order to estimate the zero at each of 
the chosen scaling temperatures. 

5. Results 

In this section we will study the behavior of the zeros as a function of the lattice size, 
both in the critical and in the spin-glass phase. Finally, we will compute the density of 
zeros and extract the rj exponent from the analysis at the critical temperature and the 
Edwards- Anderson order parameter from the scaling in the low-temperature phase. 
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Table 


2. Fits of the zeros to 


e,{L) ^ A,L- 


-^1, for L > 


As we 


can see, with 




-^min ~ 


= 8 the 


per degree 


of freedom is 


acceptable only for j — 


1,2, but with 




^min 


12 aU the 


zeros have a 


reasonable fit 


However, the value of Xi 


grows with j , 




an indication that 


we have to consider corrections to scaUng (see 


text). 




j 


-^min 




/3 








xVd.o.f. 


1 


8 




0.902 




2.703(12) 




1 'TO 1 

1.78/3 


2 


8 




0.902 




2.712(6) 




3.23/3 


3 


8 




0.902 




2.718(5) 




8.12/3 


4 


8 




0.902 




2.725(5) 




15.1/3 


1 


12 




0.902 




2.695(14) 




1.27/2 


2 


12 




0.902 




2.715(8) 




2.95/2 


3 


12 




0.902 




2.731(7) 




2.19/2 


4 


12 




0.902 




2.745(7) 




2.49/2 



5.1. Scaling at the Critical Point 

We first consider the scaling at the critical point and use it to determine the anomalous 



dimension, as in (10). Our simulations were optimized to investigate the low- 
temperature phase, for large system sizes, rather than to obtain the critical parameters. 
Therefore, we take the value of /3c = 0.902(8) from [32J, which features many more 
samples and small sizes to control scaling corrections but does not reach the low- 
temperature phase, and will also use this reference to check our value of 

Let us first consider a fit of the individual zeros, leaving aside corrections to scaling. 



i.e., following (10). For the j'-th zero, we fit to 



e^iL) = AjL-'^K (22) 
We report the results of these fits in table |2] We see that the first and second zeros 



follow (22) very well for L > 8, but for j > 2 we need to restrict the fit to L > 12. 
However, there is an inconsistency in the results: the value of Xi should be the same 
for all zeros, but we see that it increases with j. Moreover, at least for the larger j, 
Xi is incompatible with the expected value of Xi = 2.688(5), (taking rj = —0.375(10) 



from [32]) This hints that corrections to scaling should be taken into account, as in (12) 



In order to do so, we consider all values of j at the same time and perform a global 
fit, enforcing data from different zeros to share the same xi and X2- As points coming 
from a given L are correlated, the full covariance matrix has to be considered. We label 
our set of points {ej{La)} by their L and their j: we have data for £ = 5 different 
values of L (Li = 8, = 12, L3 = 16, L4 = 24, L5 = 32) and for j = 1,2,3,4. The 

II If we combine the critical exponents of |32j with the Janus simulations studied herein, we obtain a 
compatible value of f3c = 0.905(7) |33) . 



Numerical study of the overlap Lee- Yang singularities in the EA model 



9 



0.02 



0.01 




Figure 1. The four first zeros at /3 = /3c- In order to appreciate the scaling better, we 
show only the data for L > 12 and compare to equation (12), fixing 2:2 = w = 1-0(1) 
from [32] and performing a global fit for a common value of xi (see text). We obtain 
Xi = 2.67(6) [1], with a chi-square per degree of freedom of /d.o.i. = 5.88/7. 



appropriate goodness-of-fit estimator is, therefore, 

4 c 

i,j=l a,b=l 

where (T{ia){jh) is the covariance matrix of the set of zeros (which is block diagonal, since 
data for different L are uncorrelated) . 

Unfortunately, we do not have enough data to determine X2 and Xi at the same time 
(the resulting error in uj would be greater than 100%). Instead, we take X2 = uj = 1-0(1) 
from [32j and fit only for Xi and the amplitudes. The resulting fit for L > 12, shown in 
figure [T} gives 

xi = 2.67(6)[1], xVd-o.f. = 5.88/7, (24) 

where the error in square brackets accounts for the uncertainty in u. Our determination 
of Xi is now compatible with the expected value of xi = 2.688(5). Therefore, the scaling 
of the zeros is consistent at the critical point. 
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0.02 



0.01 




Figure 2. Scaling of tlie zeros at /? = 1.4, witli a best fit to (251 for L > 16. We 
obtain xi = 2.842(11), with xVd.o.f. = 7.34/7. 



0.02 



^1 0.01 - 




Figure 3. Scaling of the zeros at /3 = 1.2, with a best fit to (25 1 for L > 16. We 
obtain xi = 2.844(10), with xVd-o.f. = 2.89/7. 



5.2. Scaling in the low-temperature phase 

Now we consider the scaling of €j{L) in the low-temperature phase. This time, we 



expect, from (13), 



6,(L)^AL--S 



(25) 



with Xi = D. 
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Table 3. Scaling of the zeros in the low-temperature phase. For the two considered 
temperatures {(3 — 1.2, 1.4) we first show a fit without corrections to scaling for L > 16, 



that is ej{L) ~ AjL^'^^. As explained in Section 5.1 this is a global fit for the four 
zeros, considering their full covariance matrix. We then consider the same fit with 
corrections to scaling, trying different values for uj (see the text for more details). In 
all cases xi is smaller than the expected value xi — D = 3. 



16 


1.4 


- 


2.842(11) 


7.34/7 


12 


1.4 


1 


2.57(12) 


3.79/7 


12 


1.4 


3 


2.79(2) 


4.18/7 


12 


1.4 


0.255 


2.75(10) 


17.6/7 


12 


1.4 


0.39 


2.67(9) 


14.7/7 


12 


1.4 


0.65 


2.55(11) 


7.90/7 


12 


1.4 


0.79 


2.48(13) 


5.03/7 


16 


1.2 




2.844(10) 


2.89/7 


12 


1.2 


1 


2.82(5) 


10.5/7 


12 


1.2 


3 


2.84(2) 


6.90/7 


12 


1.2 


0.255 


2.80(10) 


12.3/7 


12 


1.2 


0.39 


2.80(12) 


11.9/7 


12 


1.2 


0.65 


2.81(10) 


11.3/7 


12 


1.2 


0.78 


2.81(8) 


11.0/7 



We have fitted the data for /3 = 1.2 (Figure |3|) and /3 = 1.4 (Figure |2|) to (|25|)|f| 
The results, in Table |3j show a value of xi incompatible with xi = D = ?>. We have 
also included corrections to scaling, using both a; = 1 (Goldstone-like correction) |M] 
and a; = 3 (Ising ordered correction) [3l], u = y = 0.255 (droplet) [351 EEl [37] , 
u = 9{0) = 0.39 (replicon and also l/z> which controls the scaling correction of qEAiL) 
[S]), u = 0.79 = 26(0) (twice the rephcon gl]) and u = 0.65 = ^(^ea) PI] but in 
neither case is the asymptotic xi = D behavior recovered (see Table [3]). In addition, we 
have forced the fits with xi = 3 and leaving free u and the statistical quality of the fits 
was bad. 

The origin of this discrepancy with the standard theory can be understood using 
Eq. |9} Notice from this equation that the scaling of the zeros depends strongly on the 
behavior of the non-connected spin glass susceptibility, so only with a divergence of this 
observable as the volume, we can recover xi = 3. However, for these two temperatures 

% The crossover length Lc which marks the change between the criticality induced by the critical point 
at Tc {L < Lc) and that of the spin glass phase {L > Lc) has been computed for different values of /3 
in reference [20]. In particular, we know that Lc{/3 — 1.2) ~ 6 and Lc{l3c = 1-4)) ~ 2.5. Hence all the 
data presented in this section belong to temperatures which lie deep into the spin glass phase. In other 
words, we can only see the critical effects induced by the spin-glass phase itself, which is critical, not 
those of the critical point at T^,. 
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(/3 = 1.2 and 1.4) this is not the case (see Fig. |4|. Notice that (g^) = x/^ has not 
reached the plateau asymptotic valu^ Hence at these temperatures the spin-glass 
susceptibility does not yet diverge as the volume. 



0.45 




0.35 



0.25 



Figure 4. x/^ ^ (9^) versus the lattice size for /3 = 1.2 and 1.4. Notice that none of 
the temperatures have reached the plateau asymptotic value. 



5. 3. Behavior of the integrated density of zeroes 

We will start our analysis of the integrated density of zeros by plotting this density at 
the critical point in FigjSj One can see that the largest lattices follow a pure power law 
as predicted by the theory. The slope, on a log-log scale, of this straight line should 
correspond with an exponent a2- Fitting only the L = 32 points we obtain 02 = 1.16(2) 



in good agreement with the theory a2 = 1.116(2) (using Eq. ^ 
To obtain this figure we have discarded in the fit the first zero]* 



and ri = -0.375(10)). 



+ In a spin glass phase, both the droplet as the RSB theory predict power law corrections on the lattice 
size, so the approach to the infinite volume values is really slow. 

* This phenomenon has been previously found in the literature. For example, the authors of |28] 
studied the anisotropic Ising model at the critical point and found a different behavior of the first zero 
in the study of the integrated density. This model shows a spreading distribution of the zeros in the 
fugacity complex plane. The authors suggest that the effect of this spreading distribution of the zeros 
is modifying the behavior of the first zero. We have not computed the complete distribution of the 
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For large L we should expect a good collapse of all points in the same power law 
curve: the non collapsing part (small L in the figure) is due to the presence of scaling 
corrections (which we also found in the previous sections). 



0.01 



0.001 



C5 



0.0001 



le-05 




0.0001 



0.001 



0.01 



Figure 5. Integrated density of zeros versus the zeros at the critical point. 02 
1.16(2). 



Now we will check the theoretical predictions for the integrated density of zeros in 
the broken phase, which predict a linear behavior in the perturbing parameter e. Notice 
that in our case the margins between the critical point and the broken phase are tight 
since in the infinite volume limit we will see a behavior e^ ii^ at the critical point which 
changes just below to e (of course, this is due to the value of the 77 exponent). 

In Figs [6] and [7] we show that the data nearly follow a linear behavior of the 
integrated density deep in the spin glass phase (more concretely at /3 = 1.2 and (5 = 1.4), 
in particular for L > 24. The non-collapsing part of the curve is produced by the 
presence of scaling corrections as at the critical point. 

However, it is easy to show that if the zeros do not follow (for the lattice sizes 
simulated), in the broken phase, a scaling as the inverse of the volume, then the 

complex zeros (only in the straight line ie), nevertheless, we know from reference |17) that the zeros 
spread in the magnetic field complex plane, so it is quite natural to assume that we will have a similar 
(spreading) spatial distribution of the zeros in e. Another possible explanation is that the behavior of 
the integrated density of zeros as j /L^ is only asymptotic. These anomalies affect only the lower order 
zeros. Notice that this phenomenon affects only to the pre-factor of the power law of the smallest zeros. 
We have seen in subsection |5. 1| that the first zero scales with the right power law. We thank R. Kenna 
for interesting comments regarding this behavior. 
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Figure 6. Integrated density of zeros versus the zeros for (3 = 1.2. 
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Figure 7. Integrated density of zeros versus the zeros at (3 = 1.4. 
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integrated density of zeros does not follow exactly a linear behavior, since 02 = D/xi 



We have discussed at the end of Sec. 5.2 that this lack of 1/V behavior is related to a 



susceptibility that is not yet diverging as the volume. 

In sec. 5.2 we have found an exponent xi = 2.842(11) for /3 = 1.2 and xi = 2.844(10) 
for (3 = 1.4, which implies that a2 = 1.056(4) and 02 = 1.055(3) for (3 = 1.2 and /3 = 1.4 
respectively. 
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Figure 8. Integrated density of the zeros, for the largest lattice L — 32 and the lowest 
temperature /3 = 1.4. Notice that we are almost, but not in, the linear regime. The 
data are well fitted with b = 1.068(10). 



In Fig. |8] we show the behavior of the integrated density of zeros computed for 
our largest lattice (L = 32) and lowest temperature, (/3 = 1.4). Notice the points are 
not lying on a straight line. A fit to aie"^ works well, with 02 consistent with the value 
computed from xi (02 = 1.068(10)). So we have obtained, numerically, 02 = 1.16(2), at 
the critical point which has changed to 1.068(10) in the broken phase||] 

In this situation, we cannot compute the order parameter directly from the linear 
behavior of the integrated density since we are not observing a fully linear behavior. 
Hence, we confront our numerical data for G{e) against the theoretical prediction for 
really small e, which is G{e) = (/5gEA/7r)e. It is interesting to note (see Refs. [27] and 
that we can recover the exact slope for a given lattice size if we substitute the value 



jj We can do the same analysis with the Xi exponent: we have obtained at the critical point Xi — 2.67(7), 
which should change in the broken phase to xi = 3, although we actually see with our numerical data 
xi = 2.842(11). 
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Figure 9. Histogram {N{e) versus e) for the 1000 first zeros computed for L = 32 and 
/3 — 1.4. Notice the lack of symmetry of the histogram and the presence of events for 
large values of the zeros. 



of the order parameter computed for this lattice size. We have followed this advice, and 
we show in Figs. 10 and 11 our data for G{e) at /3 = 1.2 and 1.4 showing L = 32 
data. In addition we have plotted the asymptotic slope using the order parameter (^ea) 
computed for L = 32 lattices for these two temperatures in Ref. [20]. Notice that we 
have a slow approach to the right slope, but also that the overall picture seems to be 
correct. 

In order to gain a better understanding of this effect, we have computed the density 
of zeros not with the average of the sample zeros but with the median of the probability 
distribution of the zerosjft] 

We show in Fig. |9] the histogram of the 1000 first zeros computed on the L = 32 
lattice at /3 = 1.4. Notice from this figure the asymmetry of the histogram and the 
presence of events at large values of the zeros (which induces a large and strongly 
fluctuating value of the mean). 



One can see in Figs. [10] and 11 that the median data produce an improved scaling, 
compared with those obtained from the mean, when comparing the data with the 

ft The probability distributions one usually finds in disordered systems present long tails due the 
presence of rare events, hence, the study of the median of this kind of distributions is also very useful 
(see for example [38l 133]). 
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Figure 10. Integrated density of tlie zeros, for tlie largest lattice L = 32 and 
temperature /3 = 1.2 using the average of zeros. We have also plotted the median 
values. We have marked the expected slope at the origin, using the Edwards- Anderson 
order parameter computed in Ref. [20] for the L = 32 lattice. 



analytical prediction (slope provided by ^ea)- 

For the sake of completeness, we can cite that the integrated density of zeros using 
the medians does not behave completely linearly but with a law e^ '^^^^) ^^qj, p _ 

6. Conclusions 

By studying the complex singularities linked with the overlap we have obtained a clear 
picture of the critical region and of the low temperature phase fully compatible with 
that obtained by other more standard approaches. 

In particular, we have studied the behavior of the individual zeros as well as the 
integrated density at the critical point. In both cases we have obtained good values for 
the T] exponent and we have seen that the data are compatible with the corrections to 
scahng pubhshed in the literature [32] ■ 

Finally, we have checked the scaling laws in the spin-glass phase, obtaining strong 
scaling corrections as found previously [20]. In addition we have obtained, by monitoring 
the behavior of the integrated density, a compatible picture using the zeros with that 
obtained from the order parameter of the model (^ea) computed in finite volumes with 
standard methods. We have also shown that the use of the median instead of the mean 
improves the overall picture. 
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Figure 11. Integrated density of the zeros, for tire largest lattice L = 32 and lowest 
temperature (3 — lA using the average of zeros. We have also plotted the median 
values. We have marked the expected slope at the origin, using the Edwards- Anderson 
order parameter computed in Ref. [20] for the L = 32 lattice. 
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